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The total cross sections for the '^®^Gd(p, 7 )'^®®Tb and ^®^Gd(p,n)'^®^Tb reactions have been 
measured by the activation method at effective center-of-mass energies 3.47 < < 7.94 MeV 

and 4.96 < < 7-94 MeV, respectively. The targets were prepared by evaporation of 30.6% 

isotopically enriched ^®^Gd oxide on aluminum backing foils, and bombarded with proton beams 
provided by a cyclotron accelerator. The cross sections were deduced from the observed 7 -ray 
activity, which was detected off-line by a HPGe detector in a low background environment. 
The results are presented and compared with predictions of statistical model calculations. This 
comparison supports a modified optical proton-|-'^®^Gd potential suggested earlier. 


PACS numbers: 25.40.Lw, 26.30.-k, 27.70.+q 

I. INTRODUCTION 

Only the lightest elements (hydrogen, helium, and 
traces of lithium, beryllium) originate from a hot Big 
Bang phase in the early stage of the Universe EH- 
All other nuclei have been produced thereafter in stars 
and stellar explosions. Most nuclides above the Fe 
group (mass numbers A > 60) originate from neutron 
captures and subsequent /3 decays and thus are called 
neutron-capture elements Two main processes 

have been identified, the s process (slow neutron capture 
process, with a main component occurring in medium 
mass Asymptotic Giant Branch (AGB) stars and a weak 
component stemming from massive stars) and the r 
process (rapid neutron capture process) [^. The 
site of the latter process is disputed but core-collapse 
supernova explosions and/or neutron star mergers may 
provide viable alternatives [T^ . There are a number of 
proton-rich isotopes between Se and Hg which cannot 
be produced by either process. They are called p nuclei 
and further nucleosynthesis processes have to be invoked 
to explain their existence [il|. Among the originally 35 
p nuclei Si, large s process contributions to ^®^Er, 
^®^Gd, and ^®°Ta were found in modern AGB models 
[ 1 ^ . Also the abundances of ^^^In and ^^®Sn may receive 
contributions from the s and/or r process [T3 |. 

The currently accepted production mechanism for the 
bulk of p nuclei is the so-called 7 process, synthesizing 
proton-rich isotopes by photodisintegration reactions 
jlll [T3| . The 7 process occurs in the outer layers 
of a massive star during its core-collapse supernova 
explosion [l5l - [r^ . As an additional site, type la 
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supernovae have been suggested, where also the required 
temperatures are achieved |18l - [2ll| . In both sites, quickly 
expanding zones of hot matter show production of 
p nuclei when the reached peak temperatures are in 
the range 2 < T < 3.5 GK. In order to produce the 
heavier p nuclei it is necessary to stay at the lower 
end of this temperature range because they would be 
completely destroyed at higher temperatures. On the 
other hand, the higher end of the temperature range is 
required for the production of the light p nuclei because 
lighter nuclei are more strongly bound and less easily 
photodisintegrated [mil 113. 

The very rare ^^®La cannot be produced in a 7 process 
but it was suggested to be formed through neutrino 
reactions on stable nuclei in core-collapse supernovae, the 
so-called iz process [H [ 23 . 

A 7 process starts with ( 7 ,n) reactions on stable 
nuclei already present in the stellar plasma. When the 
temperatures remain high for a sufficiently long time, 
proton-rich isotopes are reached for which ( 7 ,p) or ( 7 , 0 ;) 
reactions are faster than ( 7 ,n). This leads to a deflection 
of or branching in the nucleosynthesis path [il [il . Due 
to the nuclear structure influencing the binding energies, 
( 7 ,p) deflections are mainly found in nuclei with neutron 
number N < 82, whereas ( 7 , 0 ;) branchings are important 
in the region N >82 [H . 

Stellar photodisintegration rates used in reaction 
network calculations are usually derived from stellar 
capture rates by applying the principle of detailed 
balance [26i - l^ . It should be noted that this reciprocity 
only applies to stellar rates, involving thermally excited 
nuclei in entrance and exit channel, but not to 
laboratory capture cross sections and their respective 
photodisintegration cross sections, unless the cross 
sections are dominated by single transitions between 
initial state, compound state, and final state [^[^.Is^. 
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Reactions involving light nuclei often include only a few 
transitions but this is never the case in the regime of 
high level-densities encountered in nuclei participating 
in the synthesis of p nuclei. Although direct ( 7 , 11 ) 
measurements have been performed also for intermediate 
and heavy nuclei (e.g., [33-l35j), it has been shown that 
these cannot constrain stellar rates, as they only allow 
to study a single 7 transition from the ground state of 
the target nucleus. Its contribution to the stellar rate is 
less than 0 . 1 %, contrary to capture data which allow to 
constrain a much larger fraction (on the order of several 
tens of percent, depending on plasma temperature and 
the nuclear level structure) of the stellar capture rate 
(see, e.g., [Ill, [13, HB-iiH). The contribution of ( 7 ,p) 
and ( 7 ,a) ground-state transitions to the stellar rate is 
even lower due to the shift of the Gamow window to 
higher compound excitation energies (40j |. Interestingly, 
the ground-state contribution of the respective capture 
rates even is larger than those for (n, 7 ) at the same 
temperature due to Coulomb suppression of transitions 
with low relative energy [I 3 , [13,113 • This even applies to 
captures with negative reaction Q values on intermediate 
and heavy targets, allowing ground-state capture to 
always constrain a much larger fraction of the stellar rates 
than photodissociation measurements. 

Although capture measurements still cannot 
completely constrain a high-temperature stellar rate, 
either, they can be used to test the predictions of 
charged particle widths. Proton and a widths determine 
the stellar photodisintegration rates because 7 -process 
temperatures, although being high even for stellar 
plasmas, still translate to comparatively low nuclear 
interaction energies, below the Coulomb barrier when 
involving charged particles. This results in charged 
particle widths being smaller than the 7 widths and thus 
a dependence of the rates primarily on these charged 
particle widths rather than on the 7 widths [Til . [40l |. 
Experiments can probe the charged particle widths 
in capture reactions but the small cross sections at 
subCoulomb energies make it difficult to measure in the 
astrophysical energy range. Therefore, astrophysically 
relevant data are scarce. The few available data 
have shown considerable deviations from predictions, 
especially for reactions with a particles [lT|, IH, lil|. 
Low-energy proton captures have been found to be 
predicted more reliably, although small modifications 
to the optical proton-|-nucleus potentials have been 
suggested [45j |. 

Proton capture reaction cross sections for the 7 process 
have been measured before, e.g., by [13, l45l - [63| . Since 
(p, 7 ) above the (p,n) threshold are sensitive not only 
to the proton width but also to the neutron and 7 
widths, additional information is required to be able to 
distinguish the impact of the various widths (see also 
Sec. IIII B II) . Such information can be obtained through a 
simultaneous (p,n) measurement because the (p,n) cross 
section well above the threshold is mainly determined by 
the proton width. Some (p,n) measurements have been 


performed on s or r process seed nuclei that affect the 
abundances of the light p nuclei, e.g., "^^Ge ^|, ®^Se [47l| . 
®®Rb [ 13 ], and on a p nucleus ^^°Te [46j| . 

Here, we present a simultaneous measurement of (p, 7 ) 
and (p,n) cross sections of ^^^Gd. Although the 7 -process 
contribution to ^®^Gd is small and originates from 
( 7 ,n) reactions on other Gd isotopes, a determination 
of proton-induced cross sections at low energy of this 
proton-rich, heavy nucleus allows to test statistical model 
predictions well below the Goulomb barrier. For the first 
time, it was also possible to determine proton capture 
reaction cross sections below the (p,n) threshold for this 
nucleus. 

The i'^ 2 Gd(p, 7 )i 53 Tb and i52Gd(p,n)i52Tb activation 
measurements have been performed in an energy range 
between 3.47 and 7.94 MeV as a test of the statistical 
model predictions over a broader energy region. The 
details of the experiment are given in Sec. [TTl The 
experimental cross sections are compared to predictions 
in the Hauser-Feshbach statistical model in Sec. [ml A 
summary and conclusions are provided in Sec. IIYI 


II. EXPERIMENTAL PROCEDURE 

The reaction cross sections of ^^^Gd(p, 7 )^^^Tb and 
^®^Gd(p,n)^®^Tb were measured using an activation 
technique. The target preparation and the experiment 
were performed using the laboratory facilities at the 
Institute for Nuclear Research of the Hungarian Academy 
of Sciences (MTA Atomki). After the target preparation, 
the experiment has two stages, first the activation of 
the target and then the determination of the number 
of /3-unstable reaction products using an off-line gamma 
detection system. The measurement of proton induced 
reactions with the activation technique is discussed in 
detail in Refs. [laiig. 


A. Target preparation 

The 30.6% isotopically enriched ^®^Gd oxide in powder 
form was provided by the company ISOFLEX USA, 
Certificate No: 64-02-152-1453. For the target backing, 
high purity thin (2.5 pm) aluminum foils, obtained from 
Espi Metal Ltd, were used. The targets were prepared 
by vacuum evaporation. 

There are two methods available in the vacuum 
chamber to evaporate the target material; via the 
resistive heating method and the use of an electron 
gun. Because of the high melting point (ss 2350 °C) of 
gadolinium oxide Gd 2 03 , the electron gun was used for 
the evaporation. Before the evaporation, 40 mg Gd 203 
powder was pressed into a pellet with 6 mm diameter. 
This Gd 2 03 pellet was placed in a tantalum boat and 
directly heated by the electron beam. 

The target holder was made of aluminum with 9 holes 
(each hole has 1.2 cm inside diameter) and A1 backing 






3 


TABLE I: Decay parameters of the p+^®^Gd reaction products and measured photo-peak efficiencies at energies of the 7 
transitions, as used in the analysis. 


Reaction 

Product 

Half-life 

7 Energy 
(keV) 

7 Intensity 
(%) 

Detection efficiency 
(%) 

'®"Gd(p,7) 

153 Tb 

(2.34 ± 0.01) d 

212.00 

31.0 ± 0.2 

1.83 ± 0.09 

'""Gd(p,n) 

1525rpj^ 

(17.5 ± 0.1) h 

271.08 

8.6 ± 0.8 

1.53 ± 0.08 




344.28 

65.0 ± 0.1 

1.27 ± 0.06 


foils with the size of 1.5 x 1.5 cm^ were placed on these 
holes. The target holder was placed 7 cm above the 
tantalum boat for Gd 203 deposition. At the same time 9 
targets were produced with thicknesses varying between 
220 jjLg/cvo? and 310 jjLg/cw?. After the evaporation, the 
targets were fixed into target frames made of aluminum 
with 3.8 cm outer diameter, 1.2 cm inner diameter, and 
3 mm thickness. The target thickness was determined by 
weighing. The weight of the A1 foil was measured before 
and after the evaporation with a precision better than 
5 fig and the Gd 2 03 number density could be determined 
from the difference. The target with the thickness 
of 310 fig/crn^ was also examined by the Rutherford 
Back Scattering (RBS) method with a microprobe to 
investigate the target homogeneity. The RBS spectra 
were taken at the Van de Graaff accelerator of MTA 
Atomki with 2.0 MeV a particles using a 3 x 3 fim^ beam 
spot size and 100 x 100 /rm^ scanning size. The largest 
difference in the target thickness between two points on 
the target was found to be 7%. 


B. Activation 

In order to activate the ^®^Gd targets, they were 
irradiated with a proton beam in an energy range 
between 3.5 and 8.0 MeV at the cyclotron accelerator 
of MTA Atomki. The proton energy was changed by 
0.5 MeV steps in the laboratory frame. Figure [T] shows 
a schematic diagram of the irradiation chamber mounted 
at the end of the beam line. The beam current was 
measured with a current integrator using the entire 
chamber as a Faraday cup, after the last beam defining 
aperture. The beam current was kept as stable as 
possible at 2 fiA. The integrated current was recorded in 
time intervals of one minute and the changes in the beam 
current during the activation was taken into account in 
the data analysis. The activations for different proton 
energies lasted between 0.5 and 30 hours, based on the 
requirement of sufficient counting statistics. 

In order to suppress the secondary electrons from 
the targets, a bias voltage of —300 V was applied 
to the entrance of the target irradiation chamber. 
Backscattered protons were detected in order to monitor 
the target stability during each irradiation by using 
an ion implanted Si detector with a reduced entrance 
aperture of 0.5 mm diameter at an angle of 165° with 


Collimator Suppression 



FIG. 1: (Golor online) A schematic diagram of the irradiation 
chamber mounted at the end of the beam line at the cyclotron 
accelerator of MTA Atomki. 

respect to the beam direction, as shown in Fig. [T] 

The target taken from the irradiation chamber after 
the activation was transferred to the off-line gamma 
detection system to determine the number of reaction 
products of the proton induced reactions on ^^^Gd. 


C. Determination of the and activities 

In order to determine the number of Tb isotopes, the 
7 radiation following their electron-capture decays was 
counted with a 100% relative efficiency HPGe detector in 
a complete 47 r low background lead shielding. After the 
end of each activation, the ^®^Gd target was mounted in 
a holder placed 10 cm from the end of the detector cap. 
The 7 spectra were taken between 2.5 and 114 hours, 
depending on the counting statistics, and were stored 
regularly in every hour, thus different isotope decays 
could be analyzed. 

The activation technique can be used to measure 
the reaction cross sections of ^^^Gd(p, 7 )^®^Tb and 
^®^Gd(p,n)^®^Tb because /3-decay half-lives of the 
reaction products in both cases are long enough to 
determine their yields (Table HI. For ^^^Gd(p,n)^®^Tb, 
the product ^^^Tb has ground (^®^®Tb) and isomeric 
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FIG. 2: The 7 -ray spectrum recorded for a counting time of 3.64 hours after 2.08 hours of activation with a 7 MeV proton 
beam. The 7 transitions used in the measurements are indicated by arrows. 


states but disintegrates 78.8% through 

internal transitions and 21.2% by (3 decay [i^. The 
contribution of |.p total ^®^Gd(p,n)^®^Tb 

reaction cross section is negligible within uncertainties, 
as was verified in a calculation with the TALYS code 



For the analysis of the ^^^Gd(p, 7 )^^^Tb reaction data, 
the 212.00 keV 7 line was used, which has an emission 
probability of 31.0% in electron-capture of ^^^Tb 
(half-life 2.34 d). In the case of the ^®^Gd(p,n)^®^Tb 
reaction, 271.08 and 344.28 keV 7 lines with emission 
probabilities of 8 . 6 % and 65.0%, respectively, were used 
because they are associated with the decay of the 
ground state of ^®^Tb with a half-life of 17.5 h. The 
7 transitions and decay parameters used for the data 
analysis are summarized in Table ID Figure [5] shows the 
7 -ray spectrum recorded for a counting time of 3.64 h 
after a 2.08 h activation with a 7 MeV proton beam. The 
7 transitions used in the measurements, as given in Table 
m are also indicated in the figure. 

The 7 -detection system was the same as the one 
used in a previous experiment and is described 

in details in [^. The absolute photopeak efficiency 
calibration was performed for a detector-target distance 
of 27 cm, where true coincidence summing effects are 
negligible, using the calibrated sources, "^Be, ^^Na, ^"‘Mn, 
5^Go, 6°Co, s^Zn, i33Ba, 

Although the true coincidence summing effect could be 
expected to be small at a distance of more than 10 cm 
between the detector and the target, it has been taken 
into account using the following method. 


A ^® 2 Q(J target was irradiated at a proton energy of 
8.0 MeV in the laboratory frame and then the 7 spectrum 
of the reaction products was measured at both 27 and 
10 cm distances. Taking into account the time elapsed 
between the two countings, the ratio of the count rates 
taken at 10 cm to the one taken at 27 cm gives the 
ratio of the photo-peak efficiency at 10 cm to one at 
27 cm. In this way, an efficiency value at 27 cm can be 
normalized to the one at 10 cm covering the coincidence 
effect correction for the 7 transition used for the analysis. 
The photo-peak efficiencies of the 7 transitions used to 
identify the products of the investigated reactions are 
also given in Table HI 


III. RESULTS AND DISCUSSION 
A. Experimental cross sections 

The cross sections of ^^^Gd(p, 7 )^®^Tb and 
^®^Gd(p,n)^®^Tb have been measured at proton 
energies ranging from 3.5 to 8.0 MeV, in steps of 
0.5 MeV. This energy range corresponds to effective 
center-of-mass energies between 3.47 and 7.94 MeV. 
The effective center-of-mass energies are defined 

as the center-of-mass energies at which one half of the 
reaction yield for the entire target thickness is obtained 
' The experimental cross section results for the 
(p, 7 ) and (p,n) reactions are summarized in Tables 
El and lllli respectively, and are also shown in Fig. |31 
The (p,n) channel becomes dominant very fast above 
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TABLE II: 
reaction. 

Measured cross sections of the 

i®2Gd(p,7)i®®Tb 

-^beam 

pefF 

Cross section 

(MeV) 

(MeV) 

(pb) 

3.500 

3.471 ± 0.014 

7.6 ± 0.8 

4.000 

3.968 ± 0.012 

43.2 ± 4.4 

4.500 

4.466 ± 0.010 

235 ± 26 

5.000 

4.962 ± 0.011 

692 ± 65 

5.500 

5.460 ± 0.009 

1277 ± 133 

6.000 

5.956 ± 0.009 

1850 ± 178 

7.000 

6.951 ± 0.007 

2796 ± 298 

7.500 

7.447 ± 0.006 

3013 ± 356 

8.000 

7.943 ± 0.008 

3875 ± 887 


TABLE III: Measured cross sections of the 
reaction. 

^®^Gd(p,n)i®2Tb 

-^beam 

TTieff 

-^c.m. 

Cross section 

(MeV) 

(MeV) 

(/ib) 

5.000 

4.962 ± 0.011 

60 ± 6 

5.500 

5.460 ± 0.009 

875 ± 81 

6.000 

5.956 ± 0.009 

4530 ± 413 

7.000 

6.951 ± 0.007 

25915 ± 2365 

7.500 

7.447 ± 0.006 

45739 ± 4174 

8.000 

7.943 ± 0.008 

101314 ± 9241 


its threshold energy (4.8 MeV) and its cross section 
becomes much higher than that of the ^®^Gd(p, 7 )^®^Tb 
reaction (Fig. [3]). 

As mentioned in Sec. U the product of i52Gd(p,n) 
reaction has both a ground (^^^®Tb) and an isomeric 
state (i 52 mrp|ji Based on the results of a TALYS 
calculation the contribution of to the 

(p,n) reaction cross sections is less than 0.9%. This 
contribution can be neglected, especially because the 



FIG. 3: Experimental cross sections of ^®^Gd(p, 7 )^®®Tb and 
^®^Gd(p,n)^®^Tb reactions. 


short-lived isomer decays by high probability (78.8%) to 
the long-lived ground state of ^®^Tb, and is therefore 
included in the counting of the ground state decay. 
For the (p,n) reaction, hence, measurements of only 
^®^Gd(p,n)^®^®Tb have been performed. The measured 
(p,n) cross sections determined from two different 7 
transitions, as explained in Sec. Ill Cl were found to be 
statistically consistent with each other. As a result, 
the weighted averages of the ^®^Gd(p,n)^®^Tb reaction 
cross sections deduced from these two 7 transitions are 
presented. 

The uncertainty in the final results has been 
determined based on the propagation of partial errors: 
counting statistics ( 0 . 2 % to 21 %), decay parameters 
(0.4% to 9%), detection efficiency (5%), target thickness 
(^ 7%), and beam current integration (less than 3%). 
The uncertainties in the effective center-of-mass energies 
range between 0.08% and 0.4%; they were calculated with 
the SRIM code [O^ based on the proton energy loss in 
the targets. 


B. Comparison with Hauser-Feshbach statistical 
model predictions 

1. Sensitivity of reaction cross sections to modifications in 
the predicted total widths 

In order to gauge the dependence of the cross sections 
on the various ingredients of the calculation and to 
understand what properties can be constrained by a 
comparison of predictions and experiment, it is useful 
to first consider the sensitivities of the astrophysical 
reaction rates and cross sections. The sensitivity s of 
a cross section or rate C is extensively discussed in (40l| . 
It is defined as 


^new ^old 

Cold 

^^new ^^old 

Wold 


( 1 ) 


as is the standard in general sensitivity analysis. In 
the current context, the sensitivity measures by how 
much a rate or cross section C changes when one 
of the averaged total widths W appearing in the 
Hauser-Feshbach formalism is varied. For example, if the 
cross section changes by the same factor by which a width 
is varied, the sensitivity is s = 1 . 0 , if it does not change, 
s = 0. The sign of s provides information on whether the 
change in C is in the same direction as the change in the 
width (s > 0 ), i.e., increase with increasing width, or in 
opposite direction (s < 0 ), i.e., decrease with increasing 
width. 

Although the ^®^Gd(p, 7 ) rate is not directly important 
in the 7 process, its sensitivity to a variation of the total 
proton-, neutron-, and 7 width is shown in Fig.|4l Similar 
results are also obtained for most (p, 7 ) reactions in the 
p nuclei range below N = 82. Due to the Coulomb 
barrier, the proton widths are quickly becoming smaller 
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FIG. 4: Sensitivity of the astrophysical reaction rates for 
^®^Gd(p, 7 )'^®®Tb to a variation of total proton-, neutron-, 
and 7 -widths as function of plasma temperature T. The 
temperature range relevant for the nucleosynthesis of heavy p 
nuclei is marked by the shaded area. 


than the 7 widths and thus dominate the sensitivity. 
Therefore it is of astrophysical interest whether the 
present experimental data can confirm the predictions of 
proton widths. Figures [5] and | 6 ] show the sensitivities 
of the reaction cross sections of ^^^Gd(p, 7 )^®^Tb and 
^^^Gd(p,n)^®^Tb, respectively, to variations of predicted 
total widths as function of center-of-mass energy, within 
the same energy range as shown in Figs. [71 HI All cross 
sections are insensitive to the total a width as it is too 
small and therefore it is not shown in the figures. 

Below the (p,n) threshold, the (p, 7 ) cross section 
is only sensitive to the proton width, above it the 
importance of 7 and neutron width quickly increase 
with increasing energy. The situation is reversed in the 
(p,n) reaction cross section which becomes less and less 
dependent on the 7 and neutron widths with increasing 
energy. In both reactions, however, the sensitivities on 
the two widths act oppositely. This has to be taken 
into account when trying to reproduce data for both 
reactions simultaneously. Regarding the proton width, 
it dominates the cross section sensitivity at the lowest 
measured energies in the (p, 7 ) reaction and at the highest 
measured energies in the (p,n) reaction. 

It has to be noted that different nuclear properties 
enter the calculation of the widths (for details, see 
Sec.5.4 in 0 ). For the particle widths, most important 
are optical potentials required for the calculation of 
transmission coefficients and low-lying excited states. 
In the radiation width, transitions to states at higher 
excitation energy are important and therefore this 
width is sensitive to the choice of 7 -strength function 
and nuclear level density. The latter enters because it is 
used above the last included discrete excited state in a 
nucleus (see, e.g., Eq.64 in [s^). 

Glose to stability, nuclear spectroscopic information is 
abundant but nevertheless it is not trivial to decide at 



FIG. 5: Sensitivity of the cross sections of ^®^Gd(p, 7 )^®®Tb 
to a variation of total proton-, neutron-, and 7 -widths as 
function of center-of-mass energy Ec.m.- 



FIG. 6 : Same as Fig.[5]but for i®^Gd(p,n)i®2Tb. 


which excitation energy to set the cut-off for inclusion of 
experimentally determined nuclear levels in each nucleus. 
It is essential for the correct prediction of particle widths 
to use a complete level scheme. Towards higher excitation 
energies, more levels are missed in experimental studies 
and the level information given in the usual databases 
cannot be considered to be complete anymore. This 
may lead to misestimated widths in the calculation when 
contributing transitions are not included due to the 
missing excited states. In this case it is advantageous 
to only include the lowest experimental levels and use a 
theoretical level density above them, even when further 
levels have been identified at higher excitation energies 
(see, e.g., Sec. 5.4.2 in l3dl. A practical example for 
this has been discussed in |68l| . This is not an issue in 7 
widths close to stability because the mainly contributing 
7 transitions involve states at high excitation energies, in 
the unresolved resonance region, anyway [s^. 

The theoretical calculations shown in Sec. IIIIB 21 
have been obtained with version 0.8.4s of the nuclear 
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reaction code SMARAGD , which makes use of level 
schemes from [^, also included in the 2010 version 
of NuDAT [s^. At most 40 experimental states are 
used for each nucleus, as long as they are found as a 
consecutive sequence of levels with known spin and parity 
assignment. A peculiar situation arises for the neutron 
widths in the reactions studied here: the (p,n) energy 
range measured here extends into a region of incomplete 
level information because only rotational bands with 
spins J > 8 are known above 345 keV excitation energy 
in ^^^Tb. Varying neutron- and y-widths, as well 
as the included ^®^Tb levels, it was found that a 
simultaneous reproduction of the (p,n) and (p, 7 ) data 
is only possible when using a theoretical nuclear level 
density above 345 keV. This indicates that further, 
unidentified low-spin levels must be present in addition to 
the experimentally determined rotational bands, which 
is not surprising as the total level density is expected to 
strongly increase with increasing excitation energy. No 
such problem was found for the proton width, indicating 
that the relevant excited states at low energy are already 
included in the experimental level scheme. 

It should further be noted that the above treatment of 
excited states was also used to calculate the sensitivities 
shown in Figs. |H- E) Therefore they differ from the 
sensitivities given in [^, which used the default NuDAT 
set of experimental states. 


2. The optical Gd potential 

As pointed out above, the low-energy proton width is 
the quantity of astrophysical interest here. As shown 
in Figs. [5] and |6] the cross sections are sensitive to 
the proton widths across all investigated energies. The 
additional dependence on the neutron- and y-widths 
can be addressed by simultaneously comparing to the 
(p,y) and (p,n) data which exhibit different sensitivities, 
as discussed above. Since the sensitivities are almost 
symmetric in the neutron- and y widths, only the ratio 
between the two widths can be determined. For the 
reactions discussed here, however, the neutron width is 
already well determined through the requirement of (p,n) 
cross section reproduction close above the threshold, 
which was also used to study the excited states cut-off 
ISec. IIII B ip . Any remaining discrepancies between 
theoretical and experimental cross sections at the high 
end of the measured energy range of the (p,y) reaction 
and close to the (p,n) threshold are likely due to 
deficiencies in either the theoretical description of the 
proton width through an optical potential or the y width, 
which includes two only theoretically known quantities, 
the y-strength function and the nuclear level density. The 
action of the two widths can be distinguished, however, 
through their different impact on the (p,y) and (p,n) 
cross sections. 

We compared the measured cross sections to 
SMARAGD predictions using different popular optical 



Ec'm. (MeV) 


FIG. 7: Experimental cross sections (exp) of ^®^Gd(p, 7 )^®®Tb 
as function of effective center-of-mass energy , compared 
to statistical model calculations using different optical 
p-|-^®^Gd potentials: the potentials by [^. [t^I (Lej), by [t^ 
(Lane consist.), and by [33, (mod Lej). 

p-bnucleus potentials while keeping all other ingredients 
to the calculations fixed. The results are shown in Figs. [7] 
and in Essential is a simultaneous reproduction of (p,y) 
and (p,n) data. Of special interest is the comparison 
with the three lowest (p,y) data points which are below 
the (p,n) threshold, where the cross sections are only 
sensitive to the proton width. At all other energies, the 
impacts of neutron-, proton-, and y-widths have to be 
disentangled. The semi-microscopic optical potential of 
[nl , including low-energy modifications for astrophysics 
by 0 (marked ’Lej’ in Figs. [3 0), was also used in the 
large-scale reaction rate calculations of [2^. It results 
in a different energy dependence than observed in the 
low-energy (p,y) cross sections although the absolute 
magnitude is reproduced well at these energies. There is 
also a problem at the higher energies. Since the predicted 
cross sections are below the data for both reactions, it is 
not possible to amend this problem by changing the y- 
and/or neutron width. 

Similarly to the potential by |Iil[73 , the more recent 
potential by [t^ (marked ’Lane-consist.’ in the figures) is 
the Lane-consistent version of a potential obtained from 
a Briickner-Hartree-Fock model using a Local Density 
Approximation but with parameters fitted to recent 
experimental data. Although the (p,n) cross sections at 
higher energy are reproduced well with this potential, 
the low-energy (p,y) cross sections are significantly below 
the data. Ghanging the ratio of the neutron- and y-width 
cannot improve the overall agreement as it mostly affects 
the cross sections at the highest energies in the (p,y) 
reaction and the ones at the lowest energies in the (p,n) 
reaction, while the (p,y) cross sections below the (p,n) 
threshold remain unchanged. 

Finally, it was found in previous work [s^, IH, 
that an improved description of low-energy (p,y) data is 
possible with a modified imaginary part of the potential 
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FIG. 8: Same as Fig.[7|but for i®^Gd(p,n)i®2Tb. 


by [72|. Also for the two reactions discussed here, use of 
this potential (marked as ’mod Lej’ in the figures) yields 
the best overall description of the data. A further slight 
improvement at low (p,n) and high (p, 7 ) energies can 
be achieved by increasing the 7 width by about 10 % but 
this does not change the conclusions regarding the proton 
optical potential. 


IV. SUMMARY AND CONCLUSIONS 

The reaction cross sections of ^®^Gd(p, 7 )^®^Tb 
and ^®^Gd(p,n)^®^Tb have been measured by the 
activation method at effective center-of-mass energies 
3.47 < Ef^ < 7.94 MeV and 4.96 < Ef^ < 7.94 MeV, 


respectively, in order to extend the experimental 
database towards the heavier mass region for 
astrophysical reactions and to test the reliability of 
statistical model predictions. For the first time, (p, 7 ) 
cross sections below the (p,n) threshold of p+^^^Gd were 
obtained, allowing to study the prediction of proton 
widths well below the Goulomb barrier. 

Although the cross sections depend on a number of 
nuclear properties, by combining the (p, 7 ) and (p,n) 
data it was possible to disentangle the contributions of 
these ingredients and to focus on a test of the optical 
p_l_i52Gd potential. The measured cross section values 
were compared to Hauser-Feshbach statistical model 
calculations using the nuclear code SMARAGD [ 6 ^. A 
good reproduction of all data across all measured energies 
can be obtained with the recently s ugg ested potential by 
[ 45 II which is a modification of [^, l72l |. 

Further experiments to obtain (p, 7 ) data at even lower 
energies than studied here would be desirable but will 
prove very challenging due to the tiny reaction cross 
sections. 
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